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Abstract. Recent work on leading order multiparton calculations for hadronic collisions is reviewed, with 
special emphasis on the ALPGEN event generator. Some problems connected with the interface of the partonic 
events generated via matrix elements with the showering are addressed. 

PACS. 12.38.Bx Perturbative calculations - 13.85.Hd Inelastic scattering: many-particle final states 



1 Introduction 

In high energy hadronic collisions multijet final states char- 
acterize a large class of important phenomena both within 
and beyond the SM. These multiparticle final states can 
originate from hard QCD radiation processes, decay of 
SM massive particles (W, Z, top quark, Higgs boson), 
decay of heavy supersymmetric or more exotic particles. 
In general there are two different approaches to simulate 
multijet final states: the first one consists in generating 
the simplest possible final state by means of matrix el- 
ements and producing additional jets by parton shower- 
ing (HERWIG P, ISAJET 0, PYTHIA gj). This procedure 
works well in the soft/collinear regions but fails to describe 
configurations with several widely separated jets. A com- 
plementary strategy is to generate high-multiplicity par- 
tonic final states by means of exact matrix elements and 
eventually pass the generated events to further showering. 
Even if it is a leading-order (LO) approach, this proce- 
dure can become very difficult because of the complex- 
ity of the matrix element calculations with many external 
legs and of the phase-space integration. Recently there has 
been extensive activity in developing several parton-level 
Monte Carlo (MC) event generators, such as ACERMC g], 
ALPGEN 0, AMEGIC++ 0, CompHEP [7], GRACE 
HELAC/PHEGAS/ JETI [3], MADEVENT 10 . In this contribu- 
tion the state of the art of the ALPGEN generator is re- 
viewed, paying attention to the latest improvements. The 
general problem of interfacing a LO partonic event gener- 
ator with the parton shower is addressed, reviewing some 
recent work on the subject. 



2 The ALPGEN event generator 

The ALPGEN library is a collection of MC codes dedicated 
to many processes relevant to high energy hadron-hadron 

* Based on work done in collaboration with M.L. Mangano, 
M. Moretti, R. Pittau and A.D. Polosa. 



collisions. The calculations are based on partonic events 
generated by means of exact LO matrix elements, ob- 
tained with the ALPHA algorithm for assigned kinemat- 
ics, spin, flavour and colour configurations. The generated 
unweighted events, stored according to the Les Houches 
Accord #1 H21 format, are ready for the HERWIG/PYTHIA 
evolution from partons to hadrons. Up to now the avail- 
able final states in the ALPGEN package are: 

- W + N jets, Z/j* + N jets, TV < 6, 

- WQQ + N jets, Z/~/*QQ + N jets (Q = c,b,t), N < 4, 

- W + c + N jets, N < 5, 

- n W + m Z + p Higgs + N jets, n + m +p < 8, N < 3, 

- to 7 + N jets, N + to < 8 and TV < 6, 

- QQ + N jets, (Q = c, b, t), N < 6, 

- QQQ'Q' + N jets, (Q, Q' = b, t) , N < 4, 

- TV jets, N < 6, 

- QQH + N jets, (Q = b,t), N < 4. 

The limitations in the maximum number of jets is only due 
to the setting of internal parameters in the ALPHA code, 
which could be changed to accomodate a larger number of 
final-state particles. While in the first version of ALPGEN 
the top quarks were considered as real particles, they are 
now (vl.3, in the QQ and QQH processes) allowed to 
decay in the three final state fermions (t — > Wb — ► bff) 
with exact matrix element, thus retaining all the spin cor- 
relations among the top decay products. The decay is cal- 
culated in the approximation of on-shell top quark and 
W boson, in order to avoid the inclusion of non-resonant 
diagrams while preserving the gauge invariance of the cal- 
culations. The same strategy has been implemented for 
the decay of the gauge boson in the vbjets code, the 
ALPGEN directory dedicated to multiboson plus jets pro- 
duction, where the matrix clement for the vector boson 
decay into a fermionic pair has been introduced in the zero 
width approximation. The generation of multiboson final 
states requires a careful treatment of the widths in the 
propagators of the unstable particles, because they gen- 
erally break gauge invariance, giving rise to a bad high 



2 



Fulvio Piccinini: Improving Multiparton Monte Carlo Tools in Hadronic Collisions 



energy behaviour of the cross sections |13) . The strategy 
adopted in ALPGEN is to calculate the matrix element with 
the bosonic widths set to zero, removing the events con- 
taining a vector boson with a propagator mass Mq such 



that 



\M 2 - 

M 2 -s 



M 2 \ < s , with 



ds 



(s - M 2 ) 2 



Ai- 



ds 



1 



(s- M 2 ) 2 +r 2 M 2 



and 



\M -M\< 



r 



For all other processes that do not involve the presence 
of several gauge bosons, the fixed-widths prescription is 
used. The Higgs bosons are treated as stable particles. 
Their decay to fermion pairs or to four fermions will be 
implemented in a future release of the programme. 



3 Matching partonic event generators to 
parton shower 

In order to simulate the real hadronic final states, the 
partonic events have to be passed through a MC parton 
shower. However, in this interface there is some ambiguity 
in the cuts implemented at the partonic level. They are un- 
physical, in the sense that the final jet cross-section should 
be independent of their choice, provided that they are not 
harder than the cuts applied to the real jets. However, 
starting with looser partonic cuts increases the probability 
of obtaining n jets from n+m partons after parton shower- 
ing (giving rise to the double-counting problem) , as can be 
seen in fig.Q which shows the jet rates with the constraint 
E^ 1 > E 1 ™* for the hardest jet in FT + 3 jets events at 
Tevatron, versus the parton separation Z\i? part imposed at 
the level of ME generation. The jets are reconstructed with 
the cone algorithm and the cross sections are normalized 
to the result obtained with Z\i? part = 0.7. The increasing 



ratio for smaller AR part is due to both the collinear di- 
vergence of the matrix element for Z\_R part — > and the 
increasing double counting for smaller Z\i? part . A first ap- 
proximation to the solution of the problem could consist 
in requiring a jet matching for every parton |14j . With this 
recipe the shapes of the W + 3 jets example above become 
flatter, as diplayed in fig.0 but still showing a residual de- 
pendence on the parameter Z\i? part . The general problem 
of matching multiparton matrix element calculations with 
parton showers has been extensively studied in the litera- 
ture and for e + e~ collisions a solution (CKKW) has been 
proposed and tested on LEP data which avoids dou- 
ble counting and shifts the dependence on the resolution 
parameter beyond next-to-leading logarithmic (NLL) ac- 
curacy. The method consists in separating arbitrarily the 
phase-space regions covered by matrix element and par- 
ton shower, and use, for all parton multiplicities, vetoed 
parton showers together with reweighted tree-level matrix 
elements by means of suitable Sudakov form-factor com- 
binations. The necessary steps for the implementation of 
the procedure can be summarized as follows: 
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Fig. 1. The rate for pp — > VF + 3 jets at Tevatron as a function 
of the partonic separation cut zA_R pa rt normalized to the cross 
section for AR paTt — 0.7. 
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Fig. 2. The same as in fig.Q requiring a jet matching for every 
parton |14|. 



select the jet multiplicity n according to the jet rates 
obtained with matrix elements with resolution > 
Hcutt defined according to the fc^-algorithm JSj (yy = 
2mm(E 2 ,E 2 )/s(l-costi lJ )); 

generate n parton momenta according to the matrix el- 
ement with fixed a s (y cut ) and reweight the event with 
the probability of no further branching by means of 
Sudakov form factors; 

build a "PS history" by clustering the partons to de- 
termine the values at which 1,2,.. .,n jets are resolved. 
In so doing a tree of branchings is constructed and the 
nodal scales characteristic of each branching are used 
to reweight the event with running a s ; 
apply a coupling constant reweighting factor ot s {yi) / 
cts(y C ut) < 1 for every branching of the "PS history", 
where yi is the nodal scale; 
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— after successful unweighting, use the n-parton kine- 
matics as initial condition for the shower, vetoing all 
branchings such that j/y > y C ut- 

The extension of the procedure to hadronic collisions has 
been proposed in ref. Recent detailed results of the 
implementation of the procedure with the programs HERWIG. 
PYTHIA and SHERPA have been presented in ref. [H]. The 
implementation in ALPGEN is currently under investiga- 
tion. Preliminary studies regarding the partonic steps with 
ALPGEN have been presented in ref. [T^]. In the case of 
hadronic collisions there is a certain degree of arbitrari- 
ness, such as the choice of the Sudakov form factors, the 
choice of the scale of a s (LO or NLO), the treatment of 
the highest- multiplicity matrix element, the choice of the 
clustering scheme, the use of flavour or colour information 
to define the tree and the related reweighting factors, the 
treatment of gauge bosons. All these uncertainties entail 
that a large degree of tuning on the data (possibly process- 
dependent) will be needed, and further work remains to 
be done to find what the correct prescriptions are. The 
a s reweighting of the partonic events could be important 
on its own, because it should effectively give, in a gauge- 
invariant way, the bulk of the NLO QCD corrections. This 
could be tested in cases where multijet NLO calculations 
are available. 

4 Summary 

The MC simulation of hard multiparticle final states at 
hadronic colliders is a very important issue. Thanks to re- 
cent efforts by different groups, several multiparton event 
generators based on exact matrix elements are now avail- 
able; they were thoroughly cross-checked during the MC 
workshops held at FNAL and CERN during 2003. These 
programs generate samples of unweighted events in a stan- 
dardized format (the LesHouches Accord # 1) which can 
be passed to the parton shower-based MC programs to 
go from the partons to the real final-state hadrons. The 
matching between a LO multijet event generator and the 
parton shower MC suffers from the serious problems of 
double-counting and dependence on the parton-level cuts. 
For the case of e + e~ collisions the problem has been solved 
beyond NLL accuracy with the CKKW algorithm. This 
can be extended to hadronic collisions, but the proof is 
still pending. However, recently, there has been an in- 
tense activity in its implementation on existing MC event 
generators going through many subtleties involved in the 
CKKW algorithm for hadronic collisions. 

The author wishes to thank the conveners of the Hard 
QCD Parallel Session for the invitation. M.L. Mangano 
and A.D. Polosa are gratefully aknowledged for many use- 
ful discussions and a careful reading of the manuscript. 
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